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Mixing and CP-violation in charm 
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The motivation most often cited in searches for D° — D° mixing and CP-violation in charm system lies with 
the possibility of observing a signal from new physics which dominates that from the Standard Model. We review 
recent theoretical predictions and experimental constraints on D° — D° mixing parameters. We also discuss the 
current status of searches for CP-violation in charmed meson transitions, as well as some recent theoretical ideas. 



1. Introduction 

Charm transitions play a unique dual role in 
the modern investigations of flavor physics. It 
provides valuable supporting measurements for 
studies of CP-violation in i?-decays, such as form- 
factors and decays constants, as well as out- 
standing opportunities for indirect searches for 
physics beyond the Standard Model (SM). It 
must be noted that in many dynamical models 
of new physics the effects of new particles ob- 
served in s, c, and b transitions are correlated. 
Therefore, such combined studies could yield the 
most stringent constraints on their parameters. 
For example, loop-dominated processes such as 
D°—D° mixing or flavor-changing neutral current 
(FCNC) decays are influenced by the dynamical 
effects of down-type particles |Tj, whereas up-type 
particles are responsible for FCNC in the beauty 
and strange systems. Finally, from the practi- 
cal point of view, charm physics experiments pro- 
vide outstanding opportunities for studies of new 
physics because of the availability of large statis- 
tical samples of data. 

The basic idea behind searches for new physics 
in the relatively low energy processes like charm 
transitions stems from the fact that the effects of 
new interactions in can be naturally written in 
terms of a series of local operators of increasing 
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dimension. Together with the one loop Standard 
Model effects [5] they can contribute to AC = 1 
(decays) or AC = 2 (mixing) FCNC transitions. 
For example, in the case of D° — D° mixing these 
operators generate contributions to the effective 
operators that change D° state into D° state, 
leading to the mass eigenstates 

\D l2 )=p\D )±q\D ), (1) 

where the complex parameters p and q are ob- 
tained from diagonalizing the D° — D° mass ma- 
trix. Note that \p\ 2 + \q\ 2 = 1. If CP-violation 
in mixing is neglected, p becomes equal to q, 
so |.Di 2) become CP eigenstates, CP\D±) = 
±\D±). The mass and width splittings between 
these eigenstates are given by 

m 2 - mi r 2 - Ti 
x= , v= . (2) 

r ' y 2r w 

It is known experimentally that D° — D° mix- 
ing proceeds extremely slowly, which in the Stan- 
dard Model is usually attributed to the absence 
of superheavy quarks destroying GIM cancella- 
tions. This situation is an exact opposite to what 
happens in the B system, where B a — B° mix- 
ing measurements are used to constrain top quark 
couplings [3]. 

It is instructive to see how new physics can af- 
fect charm mixing. Since the lifetime difference y 
is constructed from the decays of D into physical 
states, it should be dominated by the Standard 
Model contributions, unless new physics signif- 
icantly modifies AC = 1 interactions. On the 
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contrary, the mass difference x can receive con- 
tributions from all energy scales. Thus, it is 
usually conjectured that new physics can signifi- 
cantly modify x leading to the inequality x 3> y 2 ■ 

Another possible manifestation of new physics 
interactions in the charm system is associated 
with the observation of (large) CP-violation. This 
is due to the fact that all quarks that build up the 
hadronic states in weak decays of charm mesons 
belong to the first two generations. Since 2x2 
Cabbibo quark mixing matrix is real, no CP- 
violation is possible in the dominant tree-level 
diagrams which describe the decay amplitudes. 
CP-violating amplitudes can be introduced in the 
Standard Model by including penguin or box op- 
erators induced by virtual 6-quarks. However, 
their contributions are strongly suppressed by 
the small combination of CKM matrix elements 
V c bV* h . It is thus widely believed that the ob- 
servation of (large) CP violation in charm decays 
or mixing would be an unambiguous sign for new 
physics. This fact makes charm decays a valu- 
able tool in searching for new physics, since the 
statistics available in charm physics experiment 
is usually quite large. 

As in B-physics, CP-violating contributions in 
charm can be generally classified by three dif- 
ferent categories: (I) CP violation in the decay 
amplitudes. This type of CP violation occurs 
when the absolute value of the decay amplitude 
for D to decay to a final state / (Af) is differ- 
ent from the one of corresponding CP-conjugated 
amplitude ("direct CP-violation"); (II) CP vio- 
lation in D a — D° mixing matrix. This type of 
CP violation is manifest when = \p/q\ 2 = 
(2M 12 - iT 12 )/(2M* 2 - iT* 12 ) ^ 1; and (III) CP 
violation in the interference of decays with and 
without mixing. This type of CP violation is pos- 
sible for a subset of final states to which both D° 
and D° can decay. 

For a given final state /, CP violating contri- 
butions can be summarized in the parameter 



q_Aj_ 

P A f 



R p l 



A 



A, 



(3) 



where A f and A j are the amplitudes for D° — > / 

2 This signal for new physics is lost if a relatively large y, 
of the order of a percent, is observed HI . 



and D° — > / transitions respectively and S is the 
strong phase difference between Af and Af. Here 
4> represents the convention-independent weak 
phase difference between the ratio of decay ampli- 
tudes and the mixing matrix. Since CP-violation 
in the mixing matrix is expected to be small, we 
will often expand R^ 2 — 1 ± A m 0] . 

At present, most experimental information 
about the D° — D° mixing parameters x and 
y comes from the time-dependent analyses that 
can roughly be divided into two categories. First, 
more traditional studies look at the time depen- 
dence of D — *■ / decays, where / is the final 
state that can be used to tag the flavor of the 
decayed meson. The most popular is the non- 
leptonic doubly Cabibbo suppressed (DCS) decay 
D° — > K + ir~ . Time-dependent studies allow one 
to separate the contribution of direct doubly Cab- 
bibo suppressed transition from the one involving 
mixing D° -> W ->• K + tv~, 



T[D Q ^K+ir-]=e- Tt \A K -^\ 2 

R + V~RR m {y' cos <f> - x' sin 4>)Tt 

+ ^+x 2 mr 



(4) 



where \fRe %5 parameterizes the ratio of DCS 
and Cabibbo favored (CF) decay amplitudes, and 
x' = x cos S + y sin 6, y' = y cos S — x sin 5. Since 
x and y are small, the best constraint comes from 
the linear terms in t that are also linear in x and 
y. As follows from Eq. direct extraction of 
x and y is not possible due to unknown relative 
strong phase 8 of DCS and CF amplitudes 
This phase, however, can be measured in an ex- 
periment which produces CP-correlated pairs of 
£>-mesons, such as CLEO-c Recent experi- 
mental constraints on D — D° mixing parameters 
from D° — > K + tt~ analyses are presented in Ta- 
ble (presented Belle and FOCUS results assume 
conservation of CP). 

Second, D° mixing parameters can be mea- 
sured by comparing the lifetimes extracted from 
the analysis of D decays into the CP-even and 
CP-odd final states. This study is also sensitive 
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Table 1 



Recent measurements of 


x and y in D° - 


-> K+ir- 


Experiment 


x' 2 (95% CL) 


y' (95% CL) 




(xlO- 3 ) 


(xl0~ 3 ) 


Belle (2004) 


< 0.81 


-8.2 < y> < 16 


BaBar (2003) 


< 2.2 


-56 < y' < 39 


FOCUS (2001) 


< 1.52 


-124 < y' < -5 


CLEO (2000) 


< 0.82 


-58 < y' < 10 



The experimental values are given in ref. jZj. 



to a linear function of y via 

t(D -> K-7T+) . A m 

T(D^K+K~) - 1 = yC ° S ^~ XSm ^—- (5) 

The results from various experiments can be 
found in ref. [5j, with the grand average yep = 
(0.9 ±0.4)%. 

Time-integrated studies of the semileptonic 
transitions are sensitive to the quadratic form 
x 1 + y 2 and at the moment are not competitive 
with the analyses discussed above. 

The construction of tau-charm factories CLEO- 
c and BES-HI allows for new time-independent 
methods that are sensitive to a linear function 
of y. One can use the fact that heavy meson 
pairs produced in the decays of heavy quarko- 
nium resonances have the useful property that the 
two mesons are in the CP-correlated states [0]. 
For instance, by tagging one of the mesons as a 
CP eigenstate, a lifetime difference may be deter- 
mined by measuring the leptonic branching ratio 
of the other meson. Its semileptonic width should 
be independent of the CP quantum number since 
it is flavor specific, yet its branching ratio will be 
inversely proportional to the total width of that 
meson. Since we know whether this D(k2) state 
is tagged as a (CP-eigenstate) D± from the decay 
of -D(fci) to a final state S a of definite CP-parity 
a = ±, we can easily determine y in terms of 
the semileptonic branching ratios of D±, which 
we denote B l ± . Neglecting small CP-violating ef- 
fects, 

i ( B{(v) BUpys 

V 4 \B e _(V) B l + {V) ) ' 1 ' 

A more sophisticated version of this formula as 



well as studies of feasibility of this method can 
be found in ref. 

The current experimental upper bounds on x 
and y are on the order of a few times 10 -2 , and 
are expected to improve significantly in the com- 
ing years. To regard a future discovery of nonzero 
x or y as a signal for new physics, we would need 
high confidence that the Standard Model predic- 
tions lie well below the present limits. As was re- 
cently shown [TJ]|, in the Standard Model, x and 
y are generated only at second order in SU(3)f 
breaking, 



x , y ~ sin 2 9 C x [ST/ (3) breaking]' 



(7) 



where 9c is the Cabibbo angle. Therefore, pre- 
dicting the Standard Model values of x and y de- 
pends crucially on estimating the size of SU(3)p 
breaking. Although y is expected to be deter- 
mined by the Standard Model processes, its value 
nevertheless affects significantly the sensitivity to 
new physics of experimental analyses of D mix- 
ing 

Theoretical predictions of x and y within and 
beyond the Standard Model span several orders 
of magnitude [TI\ (see Fig.^l. Roughly, there are 
two approaches, neither of which give very reli- 
able results because m c is in some sense interme- 
diate between heavy and light. The "inclusive" 
approach is based on the operator product ex- 
pansion (OPE). In the m c 3> A limit, where A 
is a scale characteristic of the strong interactions, 
AM and Ar can be expanded in terms of matrix 
elements of local operators ^2]- Such typically 
calculations yield x,y < 10 -3 . The use of the 
OPE relies on local quark-hadron duality, and on 
A/m c being small enough to allow a truncation 
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Standard Model mixing predictions 
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Figure 1. Standard Model predictions for |x| 
(open triangles) and \y\ (open squares). Horizon- 
tal line references are tabulated in 



of the series after the first few terms. The charm 
mass may not be large enough for these to be good 
approximations, especially for nonleptonic D de- 
cays. An observation of y of order 10 -2 could be 
ascribed to a breakdown of the OPE or of dual- 
ity, but such a large value of y is certainly not a 
generic prediction of OPE analyses. The "exclu- 
sive" approach sums over intermediate hadronic 
states, which may be modeled or fit to experi- 
mental data (T3|. Since there are cancellations 
between states within a given ST/ (3) multiplet, 
one needs to know the contribution of each state 
with high precision. However, the D is not light 
enough that its decays are dominated by a few 
final states. In the absence of sufficiently precise 
data on many decay rates and on strong phases, 
one is forced to use some assumptions. While 
most studies find x, y < 10~ 3 , Refs. p3] obtain x 
and y at the 1CP 2 level by arguing that SU(3)f 
violation is of order unity. It was also shown that 
phase space effects alone provide enough SU(3) p 
violation to induce x, y ~ 10~ 2 10 . Large effects 
in y appear for decays close to D threshold, where 
an analytic expansion in SU(3) f violation is no 
longer possible; a dispersion relation can then be 
used to show that x would receive contributions 



of similar order of magnitude. 

The above discussion shows that, contrary to B 
and K systems, theoretical calculations of x and 
y are quite uncertain |14|. and the values near the 
current experimental bounds cannot be ruled out. 
Therefore, it will be difficult to find a clear indi- 
cation of physics beyond the Standard Model in 
D° — D° mixing measurements alone. The only 
robust potential signal of new physics in charm 
system at this stage is observation of large viola- 
tion of CP. 

CP violation in D decays and mixing can be 
searched for by a variety of methods. For in- 
stance, time-dependent decay widths for D — ► 
Kir are sensitive to CP violation in mixing (see 
Eq.||3Jl). Provided that the x and y are compara- 
ble to experimental sensitivities, a combined anal- 
ysis of D — ► Kir and D — > KK can yield interest- 
ing constraints on CP-violating parameters 0]. 

Most of the techniques that are sensitive to CP 
violation make use of the decay asymmetry, 



Acp(f) = 



T(D -> /) - T(D -> f) 

r(-D — > /) + r(D — > /) 



(8) 



Most of the properties of Eq. © , such as depen- 
dence on the strong final state phases, are similar 
to the ones in B-physics |15| . Current experimen- 
tal bounds from various experiments, all consis- 
tent with zero within experimental uncertainties, 
can be found in |16| . 

Other interesting signals of CP-violation that 
are being discussed in connection with tau-charm 
factory measurements are the ones that are using 
quantum coherence of the initial state. An ex- 
ample of this type of signal is a decay (D°D°) — > 
/1/2 at ^(3770) with fi and fa being the different 
final CP-eigenstates of the same CP-parity. This 
type of signals are very easy to detect experimen- 
tally. The corresponding CP-violating decay rate 
for the final states fa and fa is 



F 



1 



hh 



2RI 



(2 + .x 2 - 2/ 2 )|A /l -A /2 f 



+ y 2 )|l-A /l A /2 | 2 T fl T f2 . (9) 



The result of Eq. @ represents a slight general- 
ization of the formula given in Ref. ^7]. It is clear 
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that both terms in the numerator of Eq. 10 re- 
ceive contributions from CP-violation of the type 
I and III, while the second term is also sensitive to 
CP-violation of the type II. Moreover, for a large 
set of the final states the first term would be addi- 
tionally suppressed by SU(3) f symmetry, as for 
instance, A^^ = Xkk in the SU(3)f symmetry 
limit. This expression is of the second order in 
CP-violating parameters (it is easy to see that in 
the approximation where only CP violation in the 
mixing matrix is retained, T/j/a oc |l — B%A cx 

The existing experimental constraints |18j 
demonstrate that CP-violating parameters are 
quite small in the charm sector, regardless of 
whether they are produced by the Standard 
Model mechanisms or by some new physics con- 
tributions. Since the above measurements involve 
CP-violating decay rates, these observables are of 
second order in the small CP-violating parame- 
ters, a challenging measurement. 

It is also easy to see that the rate asymmetries 
of Eq. (|SJ) require tagging of the initial state with 
the consequent reduction of the dataset. In that 
respect, it is important to maximally exploit the 
available statistics. 

It is possible to use a method that both does 
not require flavor or CP-tagging of the initial 
state and results in the observable that is first 
order in CP violating parameters ^Hj- Let's con- 
centrate on the decays of .D-mesons to final states 
that are common for D° and D°. If the initial 
state is not tagged the quantities that one can 
easily measure are the sums 



s i = r i (t) + r i (t) 



(10) 



for i = f and /. A CP-odd observable which can 
be formed out of is the asymmetry 



a, 



CP 



(/.*) 



_ N{t) 
£, + Sj D{t) ■ 



(11) 



We shall consider both time-dependent and time- 
integrated versions of the asymmetry 1|11|) . Note 
that this asymmetry does not require quantum 
coherence of the initial state and therefore is ac- 
cessible in any _D-physics experiment. It is ex- 
pected that the numerator and denominator of 



Eq. would have the form, 



N(t) 
D(t) 



s 7= e 



— T 



[A + BT + CT 2 ] 



2e~ 



\A f \ 



A- 



(12) 



where we neglected direct CP violation in D(t). 
Integrating the numerator and denominator of 
Eq. (|TT1) over time yields 



A u CP {f) = ±[A + B + 2C] 



(13) 



where D — V / °° dt D(t). 

Both time-dependent and time-integrated 
asymmetries depend on the same parameters 
A, B, and C. The result is 



A = 



B = 



-2yVR 



A 7 



A 7 



\A f ' 2 



VI 



sin <b sin S 



\A 12 



/ + 



A 7 



cos 6 cos 5 









(p*/f- 


A 7 


1 



(14) 



C = —A. 

We neglect small corrections of the order of 
0(A m x,r fx, ...) and higher. It follows that 
Eq. l|14f> receives contributions from both direct 
and indirect CP-violating amplitudes. Those con- 
tributions have different time dependence and can 
be separated either by time-dependent analysis of 
Eq. 1)1 1H or by the "designer" choice of the final 
state. Note that this asymmetry is manifestly 
first order in CP-violating parameters. 

In Eq. H14() . non-zero value of the coefficient 
A is an indication of direct CP violation. This 
term is important for singly Cabibbo suppressed 
(SCS) decays. The coefficient B gives a combi- 
nation of a contribution of CP violation in the 
interference of the decays with and without mix- 
ing (first term) and direct CP violation (second 
term). Those contributions can be separated by 
considering DCS decays, such as D — > K^'tt or 
D — > K(*>p, where direct CP violation is not ex- 
pected to enter. The coefficient C represents a 
contribution of CP-violation in the decay ampli- 
tudes after mixing. It is negligibly small in the 
SM and all models of new physics constrained by 
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the experimental data. Note that the effect of 
CP-violation in the mixing matrix on A, B, and 
C is always subleading. 

Eq. I|14fl is completely general and is true for 
both DCS and SCS transitions. Neglecting direct 
CP violation we obtain a much simpler expres- 
sion, 



A = 0, C* = 0, 

B = — 2ysin<5sin0 \f~R 



,-r ,2 




■2 


\ A f\ + 


A 7 





(15) 



For an experimentally interesting DCS decay 
D° — > K + ir~ this asymmetry is zero in the fla- 
vor SU(3)f symmetry limit, where S — |20|. 
Since SU(3)f is badly broken in D-decays, large 
values of sin<5 jS] are possible. At any rate, re- 
gardless of the theoretical estimates, this strong 
phase could be measured at CLEO-c. It is also 
easy to obtain the time-integrated asymmetry for 
Kir. Neglecting small subleading terms of C(A 4 ) 
in both numerator and denominator we obtain 



Aqp(Ktt) — — ys'mS sin^V^R. 



(16) 



It is important to note that both time-dependent 
and time-integrated asymmetries of Eqs. (|15fl and 
Hlfty are independent of predictions of hadronic 
parameters, as both 6 and R are experimen- 
tally determined quantities and could be used 
for model-independent extraction of CP-violating 
phase 4>. Assuming R ~ 0.4% and 8 ~ 40° and 
y ~ 1% one obtains \A% P (Kir)\ ~ (0.04%) sin 0. 
Thus, one possible challenge of the analysis of the 
asymmetry Eq. Ijl6|) . is that it involves a differ- 
ence of two large rates, +„.- and Y. K -^+ , which 
should be measured with the sufficient precision 
to be sensitive to A^ p , a problem tackled in de- 
terminations of tagged asymmetries in D —> Kir 
transitions. 

Alternatively, one can study SCS modes, where 
R ~ 1, so the resulting asymmetry could be 
0(1%) sin 0. However, the final states must be 
chosen such that A^ P is not trivially zero. For 
example, decays of D into the final states that 
are CP-eigenstates would result in zero asymme- 
try (as Tf = Tj for those final states) while de- 
cays to final states like K + K*~ or p + n~ would 
not. It is also likely that this asymmetry is larger 



than the estimate given above due to contribu- 
tions from direct CP-violation (see eq. I14|) . 

The final state / can also be a multiparticle 
state. In that case, more untagged CP-violating 
observables could be constructed. For example, 
untagged studies of Dalitz plot population asym- 
metries resulting from the enantiometric interme- 
diate states were proposed in j2J to study direct 
CP-violation in Bd decays. A similar study is 
possible here as well. 

As any rate asymmetry, Eq. i|ll|) requires ei- 
ther a "symmetric" production of D° and D°, a 
condition which is automatically satisfied by all 
pp and e + e~ colliders, or a correction for D° /D° 
production asymmetry. 

In summary, charm physics, and in particu- 
lar studies of CP-violation, could provide new 
and unique opportunities for indirect searches for 
new physics. Expected large statistical samples of 
charm data will allow new sensitive measurements 
of charm mixing and CP-violating parameters. 
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